A dissipative quantum Church-Turing theorem 
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We show that the time evolution of an open quantum system, described by a possibly time dependent Liouvil- 
lian, can be simulated by a unitary quantum circuit of a size scaling polynomially in the simulation time and the 
size of the system. An immediate consequence is that dissipative quantum computing is no more powerful than 
the unitary circuit model. Our result can be seen as a dissipative Church-Turing theorem, since it implies that 
under natural assumptions, such as weak coupling to an environment, the dynamics of an open quantum system 
can be simulated efficiently on a quantum computer. Formally, we introduce a Trotter decomposition for Liou- 
villian dynamics and give explicit error bounds. This constitutes a practical tool for numerical simulations, e.g., 
using matrix-product operators. We also demonstrate that most quantum states cannot be prepared efficiently. 



One of the cornerstones of theoretical computer science is 
the Church-Turing thesis In its strong formulation it 

can be captured in the following way J^HH]: "A probabilistic 
Turing machine can efficiently simulate any realistic model 
of computation." As such, it reduces any physical process - 
that can intuitively be thought of as a computational task in a 
wider sense - to what an elementary standard computer can 
do. Needless to say, in its strong formulation, the Church- 
Turing thesis is challenged by the very idea of a quantum 
computer, and hence by a fundamental physical theory that 
initially was thought to be irrelevant for studies of complex- 
ity. There are problems a quantum computer could efficiently 
solve that are believed to be intractable on any classical com- 
puter. 

In this way, it seems that the strong Church-Turing the- 
sis has to be replaced by a quantum version [2]. Colloqui- 
ally speaking, the quantum Church-Turing thesis says that any 
process that can happen in nature that one could think of as be- 
ing some sort of computation is efficiently simulatable: 
Strong quantum Church-Turing thesis. Every quantum me- 
chanical computational process can be simulated efficiently in 
the unitary circuit model of quantum computation. 

Indeed, this notion of quantum computers being devices 
that can efficiently simulate natural quantum processes, be- 
ing known under the name "quantum simulation," is the topic 
of an entire research field initiated by the work of Feynman 
J5t] . Steps towards a rigorous formulation have been taken by 
Lloyd Ja] and many others 0] . 

Quite surprisingly, a very important class of physical pro- 
cesses appears to have been omitted in the quest for finding 
a sound theory of quantum simulation, namely dissipative 
quantum processes. Such processes are particularly relevant 
since, in the end, every physical process is to some extent 
dissipative. If one aims at simulating a quantum process oc- 
curring in a lab, one cannot, however, reasonably require the 
inclusion of all modes of the environment to which the system 
is coupled into the simulation. Otherwise, one would always 
have to simulate all the modes of the environment, eventually 
of the entire universe, rendering the task of simulation obso- 
lete and futile. We argue that the most general setting in which 
one can hope for efficient simulatability is the one of Marko- 



vian dynamics |8[] with arbitrary piecewise continuous time 
dependent control In any naturally occurring process the 
Liouvillian C determining the equation of motion 



-p(t) = £ t (p(t)) 



(1) 



of the system state p is k-local. This means that the system 
is multipartite and C can be written as a sum of Liouvillians 
each acting nontrivially on at most k subsystems. In fact, all 
natural interactions are two-local in this sense. Since we are 
interested in processes which can be viewed as a computa- 
tion, we assume that the subsystems are of fixed finite dimen- 
sion. This is arguably the broadest class of natural physical 
processes that should be taken into account in a dissipative 
Church-Turing theorem and includes the Hamiltonian dynam- 
ics of closed systems as a special case. 
In this work, we show the following. 

(i) Every time evolution generated by a fc-local time depen- 
dent Liouvillian can be simulated by a unitary quantum circuit 
with resources scaling polynomially in the system size N and 
simulation time r. 

(ii) As a corollary, we obtain that the dissipative model for 
quantum computing [11] can be reduced to the circuit model 
- proving a conjecture that was still open. 

(iii) Technically, we show that the dynamics can be approx- 
imated by a Trotter decomposition, giving rise to a circuit of 
local channels, actually being reminiscent of the situation of 
unitary dynamics. In particular, in order to reach a final state 
that is only e distinguishable from the exactly time evolved 
state, it will turn out to be sufficient to apply a circuit of Km 
local quantum channels, where 



max 



2cK 2 t 2 rb 



In 2 



(2) 



is the number of time steps, K < N k is the number of local 
terms in the Liouvillian, and b and c are constants independent 
of N, t, K, and e. Some obstacles of naive attempts to simu- 
late dissipative dynamics are highlighted, and the specific role 
of the appropriate choice of norms is emphasized. 

(iv) We also show that most quantum states cannot be pre- 
pared efficiently. 
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(v) In addition, the Trotter decomposition with our rigorous 
error bound is a practical tool for the numerical simulation of 
dissipative quantum dynamics on classical computers. 

Setting. We consider general quantum systems consisting 
of N subsystems of Hilbert space dimension d. The dynamics 
is described by a quantum master equation (Q3 with a fc-local 
Liouvillian of the form 



£ = £ £ A 

AC [Af] 



(3) 



where [N] := {1,2,..., N} and £a are strictly k-local Liou- 
villians. The subscript A means that the respective operator or 
superoperator acts nontrivially only on the subsystem A and 
we call an operator or superoperator strictly k-local if it acts 
nontrivially only on at most k subsystems. Each of the Liou- 
villians £a can be written ifioll in Lindblad form Hill 



C A = -i[if A ,-] + ^D[i Ai/1 ] 



(4) 



where T>[X]{p) := 2XpX^ - {X^X, p} and may depend on 
time piecewise continuously. In particular, we do not require 
any bound on the rate at which the Liouvillians may change. 

The propagators Tc(t, s) are the family of superoperators 
defined by 



p(t)=T c (t,s)(p( S )) 



(5) 



for all t > s. They are completely positive and trace preserv- 
ing (CPT) and uniquely solve the initial value problem 



T(t,s)=C t T(t,s) , T(s,s) = id 



(6) 



The main result, which is a bound on the error of the Trotter 
decomposition, will be somewhat reminiscent of the Trotter 
formula for time dependent Hamiltonian dynamics derived in 
Ref. 1 1311 . The main challenge comes from the fact that we are 
dealing with superoperators rather than operators. The key to 
a meaningful Trotter decomposition is the choice of suitable 
norms for these superoperators. The physically motivated and 
strongest norm is the one arising from the operational distin- 
guishability of two quantum states p and a, which is given by 
the trace distance dist(p, a) := sup 0<j4<1 tr(A(p — a)). The 
trace distance coincides up to a factor of 1/2 with the distance 
induced by the Schatten 1-norm ||-||i, where the Schatten p- 
norm of a matrix A is || A\\ p := (tr(|A| p )) 1 / p . Therefore, we 
measure errors of approximations of superoperators with the 
induced operator norm, which is the so-called (1 — > l)-norm. 
In general the (p — > <j)-norm of a superoperator T € B(B(H)) 
is defined as 114 ] 



ITI 



:= sup \\T(A)\\ q . 

\\M P =i 



(7) 



The difficulty in dealing with these norms lies in the fact 
that for p < oo the p-norm does not respect fc-locality, e.g., 



\\A <E> Inxnll i = »i||-<4||i- This problem is overcome by us- 
ing the Lindblad form of the strictly fc-local Liouvillians. In 
the end, all bounds can be stated in terms of the largest oper- 
ator norm ||-X" t || of the Lindblad operators X <G £a of the 
strictly fc-local terms. The notation X € £a means that X is 
one of the operators occurring in the Lindblad representation 
(0]i of C\. From now on we assume that this largest operator 
norm a is everywhere bounded by a constant of order 1 and, 
in particular, independent of N, i.e., a £ 0(1). 

Main result. One can always approximate any dissipative 
dynamics generated by a fc-local Liouvillian acting on N sub- 
systems, even allowing for piecewise continuous time depen- 
dence, by a suitable Trotter decomposition. The error made in 
such a decomposition can be bounded rigorously. 

Theorem 1 (Trotter decomposition of Liouvillian dynam- 



£a be a k-local Liouvillian that acts 



ics). Let £ = J2ac[n] ^a 
on N subsystems with local Hilbert space dimension d. Fur- 
thermore, let the £a be piecewise continuous in time with the 
property that a = maxA max.xeC\ sup t > H^QHoc G 0(1). 
Then the error of the Trotter decomposition of a time evolution 
up to time t into m time steps is 



L r j' -1 ' 

m. ' m. ■ 



< 



cK 2 T 2 e bT / m 



3 = 1 Ac [A] 

(8) 

where c € 0(d 2k ), b € 0(d k ), and K < N k is the number of 
strictly k-local terms £a 7^ 0. This bound holds for any or- 
der in which the products over A are taken. Tc a (t^,t^^-) 
can be replaced by the propagator Tc^ (t , t ^-^-) = 
exp(r/m£ A v ) of the average Liouvillian 



fa 



m 

T 



Tj/m 



C A dt 



(9) 



r(j-l)/m 

without changing the scaling ([8]) of the error. 



All constants are calculated explicitly in the Appendix. The 
supremum in a can be replaced by suitable time averages over 
the time steps such that ||-X't|| 0O can be large for small times. 
Before we turn to the proof of this result, we discuss important 
implications. 

Implication 1 (Dissipative Church-Turing theorem). 

Time dependent Liouvillian dynamics can be simulated effi- 
ciently in the standard unitary circuit model. 



Using the Stinespring dilation 01511 . each of the Km prop 



agators Tc a (t^-,t^-) can be implemented as a unitary U A 
acting on the subsystem A and an ancilla system of size at 
most d 2k . These unitaries can be decomposed further into cir- 
cuits ll\ of at most n = 0(log Q (1/csk)) gates from a suitable 
gate set using the Solovay-Kitaev algorithm [16] with a < 4 
such that ||[/ A — £^11°° — £ sk- Note that for pure states, we 
havei||C7 \ip)(ip | W - U I ^lli < \\U-U\\oo < £sk 

and the 1-norm is nonincreasing under partial trace. The full 
error is bounded by the error from the Trotter approximation 
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dHJ plus the one arising from the Solovay-Kitaev decomposi- 
tion, in (1 — > l)-norm bounded by KmesK- 

At this point a remark on the appropriate degree of gener- 
ality of the above result is in order. The proven result applies 
to dynamics under arbitrary piecewise continuous time depen- 
dent fc-local Liouvillians. It does not include non-Markovian 
dynamics as often resulting from strong couplings. However, 
not only this result, but no dissipative Church-Turing theo- 
rem, can or should cover such a situation: Including highly 
non-Markovian dynamics would mean to also include extreme 
cases such as an evolution implementing a swap gate that 
could write the result of an incredibly complicated process 
happening in the huge environment into the system. In such an 
intertwined situation it makes only limited sense to speak of 
the time evolution of the system alone in the first place. On the 
other hand, in practical simulations of non-Markovian dynam- 
ics, where the influence of memory effects is known, pseudo- 
modes can be included 11711 . thereby rendering the above re- 
sults again applicable. 



It has been shown recently 11811 that the set of states that can 
be reached from a fixed pure reference state by fc-local, time 
dependent Hamiltonian dynamics is exponentially smaller 
than the set of all pure quantum states. In fact, a more general 
statement holds true (see the appendix): 

Implication 2 (Limitations of efficient state generation). 

Let X? be the set of states resulting from the time evolution 
of an arbitrary initial state p under all possible (time depen- 
dent) k-local Liouvillians up to some time r. For times r that 
are polynomial in the system size, the relative volume of X? 
(measured in the operational metric induced by the 1-norm) 
is exponentially small. 

Finally, TheoremQ]also provides a rigorous error bound for 
the simulation of local time dependent Liouvillian dynamics 
on a classical computer. Even though classical simulation of 
quantum mechanical time evolution is generally believed to 
be hard in time, we have the following result. 

Implication 3 (Simulation on classical computers l(l9|). 

For systems with short-range interactions and fixed time t, the 
evolution of local observables can be simulated on classical 
computers with a cost independent of the system-size and ar- 
bitrary precision, e.g., using a variant of the time-dependent 
density matrix renormalization group method. 



To approximate the evolution of local observables in the 
Heisenberg picture, one applies the Hilbert-Schmidt adjoint, 

n'= m nAc[w] T £ A ( T ^ T ^r)< of the Trotter approxima- 
tion of the propagator to the observable. Lieb-Robinson 
bounds can be used to prove that dissipative dynamics under 
short-range Liouvillians is quasi-local and that the time evo- 
lution of local observables is restricted to a causal cone lfl9tl . 
Channels of the Trotter circuit that lie outside the causal cone 
have only a negligible effect and can be removed. The error 



for the resulting approximation to the time evolved observable 
measured in the oo-norm, is system-size independent. 

This establishes a mathematically sound foundation for 
simulation techniques based on Trotter decomposition that 
have previously been used without proving that the approx- 
imation is actually possible; see, e.g., Refs. i20ll . 

Recently, CPT maps like the local channels in the Trot- 
ter decomposition (J8j have even been implemented in the lab 
I2U|. 

Proof of theorem]]] We now turn to the proof of the main 
result. First we will find (1 — > l)-norm estimates (i) for T 
and (ii) for T~ which will be used frequently. In the next step 
(iii) we derive a product formula, which we use iteratively (iv) 
to prove the Trotter decomposition. Finally, (v) we show how 
the second claim of the theorem concerning the approximation 
with the average Liouvillian can be proven. Throughout the 
proof we consider times t > s > 0. 

(i) Because any CPT map T maps density matrices to den- 



sity matrices, we have j > 1. In Ref. 11411 it is shown 

that 



\T\ 



sup 



(10) 



for any CPT map T. Any self-adjoint operator A = A + — A- 
can, by virtue of its spectral decomposition, be written as the 
difference of a positive and negative part A± > 0. Since T is 
CPT, WTiA^W, = tr(T(A±)) = \\A\\ V hence < 1, 

and finally = 1. 

(ii) For any Liouvillian K the propagator T/c(t, s) is invert- 
ible and the inverse Tj^(t,s) = (T/c(t, s)) -1 is the unique 
solution of 



-T-(t,s)fC t , T-(«,«)=id 



(11) 



From the representation of T as a reversely time-ordered ex- 
ponential, the inequality 

11^,5)11^ <exp(| HM^dr) (12) 

follows. This can be proved rigorously with the ideas from 
Ref. 12311 (see the appendix). 

For the case where K. is strictly fc-local, we use its 
Lindblad representation and the inequality ||j4pB||i < 
Plloo||p||i||£||oo to establish ||/C|| 1 ^ 1 G 0{d k ) and hence 
|| r jc(*> a )|li_n ^ e b ^- s \ with b g 0(d k ). 

(iii) In the first step we use similar techniques as the ones 
being used for the unitary case |l3(| where differences of time 
evolution operators are bounded in operator norm by commu- 
tators of Hamiltonians. Applying the fundamental theorem 
of calculus twice, one can obtain for any two Liouvillians JC 
and C 
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T/c+c(t, s) -T K (t,s)Tc(t,s) = T K {t,s)T c {t, s) J T c (r,s) J — (T K (u, s)C t Tk.{u, s)J T k (r, s)Tjc+£.{r, s) du dr 



T K {t, s)T c (t, r)T K (u, s)[K, u ,C r }T K (r, u)T K+c {r, s) du dr 



(13) 



In the next step we take the (1 — > l)-norm of this equation, 
use the triangle inequality, employ submultiplicativity of the 
norm, and use (i) and (ii) to obtain f* J" r || [K, u , £ r ] j du dr 
as an upper bound. In the case where JC and C are strictly k- 
local ||[/C u ,£ r ]|li_n £ 0(d 2k ), which follows by the same 
arguments used in (ii) to bound ||/C|| 
is only fc-local with K terms, || [JC Ul C r ] 
at most the factor K such that 



In the case where L 
, j is increased by 



\\T K+c {t,s)-T K {t, s )T c {t, s )\\ 1 ^ l e 0((t-s) 2 e b ^d 2k K) 

(14) 

(iv) The propagator can be written as 



r 



above senses. The result can be seen as a quantum Church- 
Turing theorem in the sense that under reasonable and nec- 
essary requirements any general time evolution of an open 
quantum system can be simulated efficiently on a quantum 
computer. 
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T c (r,0) = l[Tc(Tj/m,T(j 



l)/m) . 



(15) 



Using the inequality 



HTiTa-TiTaH < H^H ||T 2 - T 2 || + - T^l ||T 2 || (16) 

and Eq. (fl4l iteratively, one can establish the result as stated 
in Eq. ©. 

(v) For any strictly fc-local Liouvillian JC the propagator 
Tjc(t, s) can be approximated by the propagator of the average 
Liouvillian, 



\T/c(t, s) - cxp( / JC r dr)\ 



\b{t - sf 



(17) 



This can be shown using the techniques described above by 
lifting the proof from Ref. Il8ll to the dissipative case (see the 
appendix). A comparison of Eq. (fTTT i with Eq. < TT~4T > shows that 
the error introduced by using the average Liouvillian is small 
compared to the error introduced by the product decomposi- 
tion and does not change the scaling of the error. 

Conclusion. In this work we show that under reasonable 
assumptions the dynamics of open quantum systems can be 
simulated efficiently by a circuit of local quantum channels in 
a Trotter-like decomposition. This channel circuit can further 
be simulated by a unitary quantum circuit with polynomially 
many gates from an arbitrary universal gate set. As a corollary 
it follows that the dissipative model of quantum computation 
is no more powerful than the standard unitary circuit model. 
The result can also be employed for simulations on classical 
computers and in the physically relevant case where the Li- 
ouvillian only has short-range interactions the simulation of 
local observables can be made efficient in the system size. It 
also shows that systems considered in the context of dissipa- 
tive phase transitions ll U 12211 can be simulated in both of the 
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Appendix 

In this appendix we elaborate on some of the technical as- 
pects of our results and give explicit expressions for all in- 
volved constants. First, we give a detailed derivation of the 
error caused by the Trotter approximation for the time evo- 
lution under a time dependent fc-local Liouvillian. Along the 
way we also derive a completely general bound for the Trot- 
ter error for arbitrary (not necessarily fc-local) time dependent 
Liouvillians, which we don't need directly for the statements 
made in the paper, but which could be of interest indepen- 
dently as the bounds of our more specialized theorem is not 
optimal in certain situations. Secondly, we present a detailed 
derivation of the error that is made when the time evolution 
under the time dependent Liouvillian is replaced by that of the 
average Liouvillian on a small time step. Finally, we prove re- 
sults on the scaling behavior of e-nets used in Implication|2]to 
argue that only an exponentially small subset of states can be 
prepared with time dependent fc-local Liouvillian dynamics in 
polynomial time from a fixed reference state. Our argument 
lifts the considerations from Ref. Ill8ll to the space of density 
matrices and the physically relevant trace distance. 



Trotter approximation for time dependent Liouvillians 

We start by giving a detailed proof that for short time in- 
tervals it is possible to approximate the time evolution of a 
fc-local time dependent Liouvillian K, + C by splitting off a 
strictly fc-local part JC and performing the time evolution un- 
der C and JC sequentially. 

Theorem 2 (Product decomposition of propagators). Let C 

and JC be two time dependent Liouvillians that act on the same 
quantum system of N subsystems with local Hilbert space di- 
mension d. Furthermore, let JC be strictly k-local and let C be 
k-local consisting of K strictly k-local terms C\. For t > s 



the Trotter error is given by 

\\T K+c (t,s)-T K {t, a)T £ (t, a)^ < (t^) 2 e t(t - s W, 

(18) 

where 



b = 2(^(2 + 4<f 
c = 2a 2 + 8a 3 d k 
a 



16a 4 d 



4j2fe 



max max sup \\X V \ 

A X6CU£ As <n<i 



We will use this theorem iteratively to bound the error 
caused by decomposing the propagators of arbitrary fc-local 
Liouvillians into the propagators of the individual strictly fc- 
local terms. 

The proof of this theorem can be presented most conve- 
niently as a series of Lemmas. From the main text (point (i) 
in the proof of the Theorem, page 3) we already know that 
completely positive and trace preserving (CPT) maps are con- 
tractive: 

Lemma 3 (Contraction property of the propagator). 

Let T be a CPT map. Then \\T\\ 1 ^ 1 = 1. 

We also need to bound the norm of the inverse propagator. 
Lemma 4 (Backward time evolution). For t > s: 

(i) Tc(t, s) is invertible and its inverse is T^(t, s) as de- 
fined by Eq. (II 11 1 in the main text. 

(ii) If the Liouvillian C is piecewise continuous in time then 



|T£(f J s)|| 1 _> 1 <exp(/ H-CrlUi dr). 

J S 



(19) 



Proof. First, we consider the case where C is continuous in 
time and use the theory presented in Ref. 112 311 and in partic- 
ular the "properties" which are proven in this reference. The 
product integral of C is defined analogously to the Riemann 
integral, 



nexp(£ r dr) := lim Texp(£ r , A 



(20) 



3=1 



where X, := XjXj- X ...X x . Since Tc(t,s) solves 

the initial value problem in Eq. © from the main text, 
Tc{t, s) = Jl! exp(£ r dr) which is exactly the statement of 
property 1. 

(i) Property 3 precisely states that a product integral is in- 
vertible. It is not hard to see that the inverse of Tc(t, s) solves 
the initial value problem (fTTl i from the main text. 

(ii) The inverse propagator is 



27(t,s) = rQexp(£ r dr) 



(21) 
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Since matrix inversion is continuous, 

1 



The definition of the Riemann integral finishes the proof for 
the continuous case. 

If £ is only piecewise continuous in time then (i) and (ii) 
hold for all the intervals where £ is continuous and from that 
and the composition property Tc{u, v)Tc(v, w) = Tc(u, w) 
We call this the reversely ordered product integral and use the (u>v>w)it follows that (i) and (ii) hold on the whole time 
convention ]T]=j Xj := X±X 2 . . . Xj. Using the submulti- interval [s, t]. □ 



T7(t,s)= lim Texp(-£ r ,Arj). (22) 

Arj— fO Vj 



plicativity of the (1 — > l)-norm and the triangle inequality we 
obtain from Eq. 



TcMW^ < lim v .n 6^(11^11^ Ar,-) (23) 



With these tools at hand we can now prove a bound on the 
Trotter error of two arbitrary (not necessarily fc-local) time 
dependent Liouvillians. 



:exp( lim > ||£ r . L . Ar 

3=1 



(24) 



J 



Theorem 5 (General Trotter error). For two arbitrary time dependent Liouvillians K, and £ the Trotter error is given by 

\\T K+C {t, s)-T K (t, 5)1^,5)111^ <J s j m^Cr)]]^ du dr e 2 O^^i d " 

<\{t-s) 2 sup || [JC U , £ r ]|li_>.i cxp(2(< - s) sup ||/C t ,|| 1 ^ 1 



t>r>u>s 



t>v>s 



(25) 
(26) 



Proof. We use a similar argument as in Ref. [13]. With the fundamental theorem of calculus we obtain 

T^(t,s)Tj^(t,s)T^ +c (t,s) -id = f d r {T c {r 1 s)T K {r,s)T K+c {r 1 s)) dr 

T c ( r ' s )i T K ( r > s )> £r]7)c+£(r, s) dr 

T £ ( r > s )( T k s)C r T K (r, s) - C r )T K (r, s)T K+c (r, s) dr 
r d 



T c( r ' S ) / ^( T K( U ' S ) £ '' T 'c( U ' S )) dM ' r K;( r '' S ) T 'C+^( r ' S ) 



dr 



T £ ( r > s ) T /c ( u i s )[£r,fc u ]Tic(u, s)T K (r,s)T K+ c(r, s) dudr 



Multiplying with T/c(t, s)T c (t, s) from the left yields 



T/c+c(t, s) - T K (t, s)T c (t, s) 



T K (t, s)T c (t, r)T K (u, s)[C r , K, U ]T K (r, u)T K+c {r, s) du dr. 



(27) 



With submultiplicativity of the (1 — > l)-norm and the bounds on the norms of the forward and backward propagators from 
Lemma [3] and [4] the result follows. □ 



I 

To complete the proof of Theorem [2] one needs to bound k-local consisting ofK strictly k-local terms £a- Then 
the norms || [C r , JC U ] \\ and ||/C r -|| 1 _ >1 in d26T i for the special 
case that K, is strictly fc-local and £ is fc-local with K strictly 
fc-local terms. 



and 



2||/C„|| 1 ^ 1 <6„ (28) 
i || [£ r ,/C„] || < c^lSr , (29) 



whereby = 4a v + 8d k a 2 , c r , u = 2a r a u +A(a r a 2 l +a 2 a u )d k + 
\%a 2 r a 2 u d 2k , and a t = maxA max{||X t || oo :XeKU £a}. 



Lemma 6. Let JC and £ be two Liouvillians which act on the 
same operator space of N subsystems with local Hilbert space 
dimension d. Furthermore, let K, be strictly k-local and £ be Proof. First, let both Liouvillians be strictly fc-local. Hence 
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each of them can be written with at most d k Lindblad opera- 
tors. Let the Lindblad representations of JC and C be 



K = -i[G,-] + J2v[K v ] 



(30) 



and 



C = -i[H r ] + ^V[L li ], 



(31) 



where V[X](p) := 2XpX^ - {X+X, p). Inequality {28]) fol- 
lows from counting the number of terms in (l30t and using that 
\\ApB\\ 1 < PH^ ilpll! llBHoo- Similarly, by writing out the 
commutator [K,,£] and using the above representations one 
can verify that [hZ ri C u ] < 2a r a u + A(a r a 2 u + a 2 a u )d k + 
ld>a 2 a 2 d 2k . If C = J2ac\n] is fc-local with K terms the 

□ 



*Ac[N] ■ 

bound is increased by at most a factor of K. 



Theorem [2] follows as a corollary of Theorem |5] and 
Lemma|6]by inserting the suprema of the bounds (|28) and ((29) 
into Eq. d26l i. Instead of using suprema in the step from 
Eq. (l25T l to Eq. d26b one can take averages over b v and cv )U 
to obtain a better, but more complicated bound. One can also 
improve the scaling of the error with the size of the time steps 
by using higher order Trotter schemes as in Ref. fl 1 311 (time 
dependent case) or Ref. 12411 (time constant case). 



Approximation by the average Liouvillian 

In the product formula in our Theorem Q] in the main text 
one can replace the time ordered integrals Tc A (t, s) by or- 
dinary exponentials of the time averaged Liouvillians. This is 
not essential to our argument concerning the quantum Church- 
Turing thesis, but makes the result more useful for applica- 
tions. The additional error caused by doing this is bounded in 
the following theorem: 

Theorem 7 (Approximation by the average Liouvillian). 

Let K, be a strictly k-local Liouvillian acting on an opera- 
tor space with local Hilbert space dimension d. Then for any 

t > s 



\\T K (t, s) - exp((t - s)/C av )lli^i < IK* - s) 2 



(32) 



where the average Liouvillian 



IC r dr 



(33) 



is indeed a Liouvillian, b = 2a 2 (2 + Ad k ), and a 
max XtG K;sup t \\X t \\ 



Proof. We lift the proof from Ref. [18] to the dissipative set- 
ting. Let t > s be fixed. Applying the fundamental theorem 
of calculus and the definition of /C av , we obtain 



T K ~(t,s)-T K (t,a) = -T K (t,8) J T^(u,s)(IC u -K. llv )T^(u, s )du 

Tfc(t, u) (JC U - K r ) T K *v (u, s) dr du 
(T)c (t, u)JC u T K ^ (u, s) — T)c(t, r)JC u T K ^ {r,sfj dr du 



1 

t-s 
1 



s J s 
t rt 



t- s 



The inequality in Eq. < fT~6b from the main text yields 



llT^ftaJ-TKMII^ < 



t - s 



ll^ulli-n (\\ T ic(t,u) -T K (t,r)|| 1 ^ 1 ||T/ C a.v(u,s)|| 1 ^. 1 
+ ll 7 *:^ Oil i^i \\Tk™(u,s) - T K av(r, s) || ^jdrdu 



(34) 



From T K {u, s) - T K {r, s) = - £* T K {v, s)KZ v dv, LemmaEl 
and the submultiplicativity of the norm we know that for 

t > u.r > s 



\T K (u,s) -T K (r,s)\\^ 



< 



l/CJ 



dv 



(35) 



and similarly for /C av . With 

Hljc^Ct, a)-2]c(*,*)||i 
< 2 



we obtain 



f f f\n^dv 

J s J s J u 



dr du . (36) 



It remains to show that /C av is a Liouvillian, i.e., that 
cxp(i/C av ) is a CPT map for all t > 0. First of all, finite sums 
of Liouvillians are Liouvillians. Furthermore, limits of se- 
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quences of Liouvillians are Liouvillians since the exponential 
function is continuous and the set of CPT maps is closed. □ 



Efficiently preparable states constitute an exponentially small 
subset of state space 

In the following we will argue that for every fixed initial 
state, the time evolution for a time interval of length r un- 
der any (possibly time dependent) fc-local Liouvillian yields 
a state that lies inside of one of Nt e-balls in trace distance. 
For times r which are polynomial in the system size, Nt is ex- 
ponentially smaller than the cardinality of any e-net (in trace 
distance) that covers the state space S. The case of Hamilto- 
nian dynamics and state vectors is investigated in Ref. Il8ll . It 
will be convenient to use the Bachmann-Landau symbols O 
and fi for asymptotic upper and lower bounds up to constant 
factors. 

By using Theorem[T]of the main text, which provides an er- 
ror bound for the Trotter approximation of a Liouvillian time 
evolution, together with the Stinespring dilation [1 1 311 and the 
Solovay-Kitaev algorithm fl6ll . one obtains the following: 



Theorem 8 (Number of channel circuits). The propagator 
from time to time t, generated by any k-local time dependent 
Liouvillian acting on N subsystems with local Hilbert space 
dimension d £ 0(1) can be approximated in (1 — > l)-norm 
to accuracy e > with one out of Nt channel circuits, where 



log(A T ) £ O 



N 3k+2 T A 



(37) 



Proof. According to Theorem[TJof the main text, the propaga- 
tor Tc(t, 0) of the Liouvillian time evolution can be approxi- 
mated by a circuit YYjLi ilAcfiV] -^A °f at most N k m strictly 
fc-local channels Tl to precision ei in (1 — > l)-norm, where 
according to Eq. (O from the main text, m = 2cN 2k r 2 / e\. 
We have assumed that 2 hx{2)cN r/e\ > b where c and b 
are given explicitly in Theorem [2] and depend only on strictly 
local properties of the Liouvillian. Employing the Stine- 
spring dilation Il5ll for each of the channels Ti one obtains 
a circuit of at most N k m strictly 3fc-local unitary gates U^. 
Each U]^ acts on an enlarged system composed of the d - 
dimensional original subsystem and an ancilla system of di- 
mension d 2k . One can use the Solovay-Kitaev algorithm lfl6ll 
to approximate every single gate XJ\ of the unitary circuit 
by a circuit \J\ of one- and two-qubit gates from a univer- 
sal gate set of cardinality 71sk £ 0(1), e.g., «sk = 3. With 
Ask = csk log Q (l/esK) of those ?t,sk standard gates, each 
unitary Ul can be approximated to accuracy esK introducing 
a total error e 2 = N k mesK- The constant c$k depends on d 3k . 

Consequently, we have for the dilation U of 
njli IIac^] -^a an approximation U with opera- 
tor norm accuracy e 2 , given by a unitary circuit of 
Aah gates = AsKA fc m standard gates from the univer- 
sal gate set. Note that for any pure state \tp), we have 



\\\U \i>){i>\U\U \i>){i>\U% < \\U - UWoo and the 
1-norm is non-increasing under partial trace. Tracing out the 
ancillas, we obtain an approximation T of Tc(t, 0) with error 
\\Tc(t,0) - < e = ei + 2e 2 . The total number of 

different channels T, which can arise in this way from the 
chosen universal gate set, is Nt < nsK NA " s ' , "" , , i.e., for given 
c,r,k,N and d, a number of Nt standard gates are enough 
to approximate any Tc(t, 0) in (1 — > l)-norm to accuracy e. 
To conclude, we bound the order of Nt- 

\og(N T ) < N m siates lognsK 

, a {2cN 3k T 2 \2cN 3k T 2 , 

= C S K log log 7JSK 



£l£2 



El 



, onai a ( 2cNt\ 2cN 3k T 2 
< c SK (3k) a log" lognsK • 

(38) 

Since we are interested in the scaling of log(Ay) for large 
N and small e\, e 2 we can assume that the argument of the 
logarithm is larger than 18 and use that log 2 (a:) < x 2 for 
x > 18 to obtain 



N 3k+2 4 

log(A T ) < C 3 2 



with C = c SK (3k) a {2c) 3 log n SK 



(39) 
□ 



The above theorem shows that the time evolution under a fc- 
local Liouvillian can be approximated by one out of Nt many 
circuits to accuracy e. The states that can be reached by any 
fc-local Liouvillian time evolution, starting from a fixed initial 
state, are hence all contained in the union of Nt e-balls (in 
1-norm) around the output states of these circuits. 

Let us now determine whether those e-balls can possibly 
cover the whole state space. For this purpose we introduce 
e-nets. We consider a D-dimensional Hilbert space % and 
denote 

(i) the set of state vectors, i.e., the set of normalized vectors 
in H by P c H, 

(ii) the set of density matrices by S C B(T~L), and 

(iii) the set of rank one projectors by V C S. 

For an arbitrary subset R C B(H) and some e > we call a 
finite subset JV[(R) C R satisfying 

WaeR3beMf(R):\\a-b\\ p <e (40) 

an e-net for R in (Schatten) p-norm. Furthermore, we call an 
e-net JVf (R) optimal if any other set X C R with smaller car- 
dinality \X\ < \J\fP(R)\ cannot be an e-net for R in p-norm. 
Similarly, we define e-nets A/" C HS (P) C P for state vectors in 
Hilbert space norm and, as before, we denote optimal e-nets 
byAC HS (P_K 

In Ref. 112.511 it was shown that for the set of state vectors of 
a D-dimensional quantum system there exist e-nets of cardi- 
nality at most |A/J iS (P)J< (5/(2e)) 2D . As the Hilbert space 
distance upper bounds 12611 the trace distance, 



\i>)-\4>)\2>k\\mii>\-\<t>)(<t>\\\i 
=diflt(ii&>w,MM), 



(41) 
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this also implies the existence of e-nets for V in p-norm of 
cardinality |JVf < (5/e) 2D for any p > 1. By comparing 
the volume of the e-balls with the volume of the whole set of 
state vectors one can see that for state vectors this construction 
is essentially optimal. 



Lemma 9. For a D -dimensional quantum system 



i\2D-l f , N 2D 

W™(P)\en((-\ )no( - 



(42) 



Proof. The set of state vectors in a D-dimensional Hilbert 
space is isomorphic to a (2D — 1) -sphere with radius 1 in 
(2£>)-dimensional real Euclidean space such that the Hilbert 
space norm | ■ I2 on state vectors coincides with the Euclidean 
norm in R 2D . The surface area of a (n — 1) -sphere of radius 
r is S n -i(r) = nC n r n -\ where C n = ir n / 2 /T(n/2 + 1) 
and r is the Euler gamma function. The set of states within 
Hilbert space distance e to a given state is a spherical cap on 
that sphere with opening angle 4 arcsin(e/2). For e <C 1, the 
area of such a cap is approximately equal to the volume of a 
(2D — l)-ball of radius e. In fact, a more detailed analysis 
reveals that for D = 3 the two are exactly identical and for 
D > 3 the cap is always smaller than the (2D — l)-ball. The 
volume of an n-ball of radius r is V n (r) = C n r n . Thus for 
D > 3, 



2D 



2DC- 



2D 



2e/ K2u-i(e) C 2 d-i£ zu 1 



^r(D + l/2) fl\ 2D ^ 15tt/1 x2 

= 2V7T — I — I 



T(D) 



where the first inequality follows from Ref. H25H . 



□ 



This is essentially the argument used in Ref. II18I1 to estab- 
lish that Hilbert space is a "convenient illusion". However, 
the lower bound on |A/" e HS (P) | does not immediately imply a 
lower bound on |A/'/'(7 : ')| (and hence also not for \JVf(S)\) 
for any p > 1. In particular, there are states with distance 2 
in Hilbert space norm and distance in any of the p-norms, 
namely, any pair of state vectors { \ip ) , — \ip )}. 

We now show that a similar lower bound as in the last 
lemma holds for the size of optimal e-nets for V and S in 
p-norm. 

Lemma 10. For p£ {1,2} 

m(S)\>\Afl(V)\en(^ 2D \ (43) 

Proof. For a given state vector \ip ) it will be convenient to use 
the notation %p := \ip)(ip\. 

We start to prove the first inequality. Fix p 6 {1, 2}. There 
is a family {pj} C J\f[(S) such that their e-neighborhoods in 
p-norm cover V and such that for each pj there exists a rank-1 



projector ipj £ V satisfying \\pj — ipj\\ p < £■ Then {i/jj} is a 
(2e)-netfor"Pinp-normwith 1^(5)1 > |{^}| > |AT e (P)|. 

From ll-H 1 > ||-|| 2 it follows that \N 2 (V)\ < |A^(P)|. 
Hence it remains to prove the lower bound for \Af 2 e (*P) I m 
j43l ). For this we construct an e'-net J\f^ s (P) for state vectors 
in Hilbert space norm from a (2e)-net M 2 t (V). For every ele- 
ment ipj £ M 2 J(P) we fix an eingenvalue-1 eigenvector \ijjj ). 
Using the (e 2 /2)-net AAi /2 ([0, 1[) = {e 2 , 2e 2 , . . . , [l/e 2 ^ 2 } 
for [0, 1[ with cyclic boundary conditions we define the set 



A/;v b (p) = { c 



27ri(5 



|V>:«e^ /2 ([0,lD, |^)e{|^)}}. 

(44) 



This is an e'-net for P and we will find an expression for e' in 
terms of e. 

Let \4>) £ P. Then there exists a state vector \tp) £ { \ tpj)} 
such that 

(2e) 2 >!|0-V'||^ = 2-2|(0|^)| 2 
>2-2|(^)|, 

and a 5 £ A/J , 2 ([0, 1[) such that 



||(^)|-Re(e 2 ^ < (2e) 2 . 

Together this yields 

3(2e) 2 > 2 - 2Rc(c 2 ^ (</>\i>)) = | |^> - e w |^) 



Since e 2,ri<5 \tp) £ Af^ s (P), we can choose e' = 4e > VT2e 
to make Af^ s (P) a (4e)-net. From the definition (04li of 
Af^ s (P) we can bound its cardinality 



1^(^)1 = l^([o,iDII{l^>}l 

< \l/e 2 ]\Ml(V)\, 



(45) 



where we have used that by construction ({IV'j)}! = 
|A/ 2 2 e ('P)|- Finally, as the described construction works for 
any (2e)-net A/ 2 2 e ('P), we obtain 



[l/e 2 ] 1^(^)1 > \Af* S (P)\ 
and Lemma|9] finishes the proof. 



(46) 
□ 



Combining Theorem|8]and LemmafTol we arrive at the fol- 
lowing theorem: 

Theorem 11 (Limitations of efficient state generation). For 

every fixed initial state, the time evolution for a time in- 
terval of length t under any k-local Liouvillian acting on 
N subsystems with local Hilbert space dimension d yields 
a state that lies inside one of Nt e-balls in 1-norm with 
log(TVr) £ O (N 3k+2 r 4 /e 5 ). For times r polynomial in the 
system size N, this is asymptotically exponentially smaller 
than \og\N}(S)\ £ n(d N /log(l/e)) where \Af*{S)\ is the 
cardinality of an optimal e-net in 1-norm that covers the state 
space S. 



